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Abstract 

The Hubbard model is reformulated in terms of different "colored" 
fermion species for the electrons or holes at different lattice sites. An- 
tiferromagnetic ordering or d-wave superconductivity can then be de- 
scribed in terms of translationally invariant expectation values for 
colored composite scalar fields. A suitable mean field approximation 
for the two dimensional colored Hubbard model shows indeed phases 
with antiferromagnetic ordering or d-wave superconductivity at low 
temperature. At low enough temperature the transition to the antifer- 
romagnetic phase is of first order. The present formulation also allows 
an easy extension to more complicated microscopic interactions. 
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The Hubbard model |L) is one of the most studied models for electron 
systems. In particular, the two dimensional model appears to be a good 
candidate @ for an explanation of high T c superconductivity. Despite its 
simplicity, several obstacles render even its approximate solution a difficult 
theoretical task. As a fermionic system it is not easily accessible to numer- 
ical simulations. Furthermore, there seems to be a competition between an- 
tiferromagnetic order and d-wave superconductivity. Both operators do not 
correspond to fermion bilinears on the same lattice site. 

Recently, promising investigations using approximations to exact renor- 
malization group equations || have been started ||, ||. The difficulty of 
these approaches, however, consists in the high complexity of the equations 
if the full momentum dependence of correlation functions for several fermions 
is kept. In particular, the low temperature phases can only be realistically 
described if effective interactions involving more than four fermions are in- 
cluded. In our opinion a prerequisite for a successful use of these methods 
in the ordered phases is a simplification of the momentum dependence of 
the interactions. This can be done if the most prominent physical degrees 
of freedom are identified. We propose here a version of the Hubbard model 
where the relevant order parameters correspond to translationally invariant 
vacuum expectation values of scalar fields. We will see that this formulation 
can describe the low temperature phases in a very simple way. We therefore 
hope that it constitutes a good starting point for a detailed renormalization 
group analysis. 

We consider the partition function of the Hubbard model |J 
Z = J D^D4>*exp{-S$,r}+ drY^ivWi + Viti)-^} (1) 

Jo L % v 

-~^*f^(^)) + (2) 

as a functional of sources r], if of the fermions as well as sources of fermion 
bilinears (Sj) that will be specified below (cf. eq. @). The spinors ipi = 

(^il, V>ij.) T , V>* = ( 1 0*t)V'ii) ( as weu as the fermionic sources rji, rj*) are two- 
component Grassmann variables depending on the Euclidean "time" r with 
antiperiodicity ipi((3) = —ipi(0), i>*(@) = — V>*(0). Here/? = 1/T is the inverse 
temperature. We treat ijj and ip* as independent Grassmann variables, even 
though the notation is reminiscent of a type of complex conjugation which 
also inverts r and reorders all Grassmann variables. (In quantum field theory 
the invariance of the action under this discrete transformation is related to 
O st erwalder- S chr ader posit ivity. ) 
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The index i labels the sites of the lattice. We concentrate on a quadratic 
lattice in two dimensions, % = (m,n), m, n G Z, with next neighbor inter- 
actions where %j = —t for i,j next neighbors and %j = otherwise. They 
describe the probability of fermion tunneling between different lattice sites. 
After a Fourier transform the Fermi surface for U — > is given by 

— 2t (cos(agi) + cos(ag 2 )) = // 

1^1 < 2A, A = n /2a (3) 

with the lattice spacing a. The localQ Coulomb interaction of the fermions 
involves the Pauli matrices f and can be rearranged, e.g. (%l)*fipi)(ip*fipi) = 
-3(ij*tpi)(tp*iJi) = -bn^nii with n iUl) = ^* Ta) ^t(l)- As usual, the expec- 
tation values of operators are related to appropriate derivatives of Z with 
respect to the sources. In addition to the discrete lattice symmetries, the 
model has two obvious continuous symmetries: the S77(2)-spin rotations and 
the [/(l)-phase rotations corresponding to charge conservation. 

In units where h — c = ks = 1, the parameters T, /i, t and U have 
dimension of mass whereas are dimensionless. The partition function 

is invariant under the rescaling (a G R + ) r — > r/a, [i — > a/j,,t at,U — > 
all, T —>■ olT . It can therefore only depend^ on the dimensionless parame- 
ter ratios fi/U,T/U,t/U and is independent of the dimensionless combina- 
tion all. Furthermore, the invariance of Z under the discrete transformation 
V>(r) — ► — — 7"), , 0*(r) — > — r ), A* ~ > ~~ ^5 t ^ —t (with appropri- 
ate transformations of the sources) permits a restriction to /x > 0. Finally, 
we may divide the lattice sites % into two classes, z G ii, if m and n are 
both even or both odd, i G J2 otherwise. The transformation ^ G / 2 — > —ipi & i 2 , 
V'ig/a "~ ~i J iei 2 while leaving ^jg/j , t/)* 6/i invariant maps Z(t) — >■ Z(—t) (again 
with an appropriate mapping for the sources). We restrict our discussion to 
positive t and //. Predictions for models with negative t or negative \x can 
easily be obtained from our results by an appropriate mapping. 

In order to represent the fermion bilinears of interest in a simple form, 
we choose to label the variables at four neighboring lattice sites by different 
colors red, yellow, green and blue, ip a , a = 1,...,4. (This also allows us 
to easily extend the formalism to lattices with different atoms.) We take 
x = (x, y) = (ma, na) with m, n even as the sites of a coarse lattice with 
lattice distance 2a and define 

lf>A(x) = ?pm,n-l , = VWl,„-l (4) 

and similar for tp* a (see fig.[j]). 

4 Here the term local refers to our neglect of the interaction of fermions located at 
different lattice sites. 

5 This holds up to a possible temperature-dependent factor from the functional measure. 
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Figure 1: The colored Hubbard model. The lattice sites of the coarse lattice 
are symbolized by a ♦. The numbering of the sites of our original lattice is 
also shown. 



The lattice symmetries T x (translation in x by a) , T y (translation in y by 
a), R (clockwise rotation by 90° around x = 0), I x (reflection at x-axis) and 
I x (reflection at the axis y = |a) act as 

T 

J. , ■ 



T 



R 



ipi(x,y 

ip<i(x,y 
ij>i(x,y 



fofa y) , ifa(x, y) -^ipx(x + 2a, y) 
">fo(x, y) , tp 3 (x, y) -> tp^x + 2a, y) 
i) 4 {x, y + 2a) , tp 2 (x, y) -> ip 3 (x, y + 2a) 
i>%{x,y) , $4(x,y) -> ipi(x,y) 

-x) , ^(x, y) -> Tp 3 (y, -x) 
^(y,-x) , ipi(x,y) ^ipi(y, -x) 
^(x,-y) , ip 2 (x,y) ^ip 3 (x,-y) 



Ix : i>i, 2 (x,y) -> i>i,2(x, ~y) , ^3,4(ar, y) -> ip 3A (x, -y + 2a). (5) 

The lattice symmetries of the coarse lattice can be composed from the gen- 
erators T x , R, I x . We note that they also act in color spaceR. 



^The local interaction is also invariant under relabelings at fixed x 



i>i -* ^2 -> if>3 -> ip4 

$1 f-» 4>4 , V>2 <-> V>3 
$1 -03 , "02 <"> V>4- 



This, however, is not a symmetry of the next neighbor interaction. 
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Useful fermion bilinear operators are (c = ir 2 ) 

a ab (x) = tpl(x)tp a (x) 
¥ab(x) = ^* b {x) r ip a {x) 

Xab (X) = $ (x) C i> a (x) = XbJ (X) 

= -toc^* r (5)=x£ ) (5). 

(We have omitted for simplicity fermion-fermion pairs in the triplet of the 
spin group, similar to the fermion- ant ifermion pairs (p. They are antisymmet- 
ric in the color indices.) Among the electron-electron (or hole-hole) pairs we 
concentrate on 



(6) 
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+ X22 + X33 


+ X44 
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-X44 
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. ~(a) . . ~(a) 
+ X23 + X34 + X41 




~(a) 
— X 2 3 


X41 > 
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= Xi 2 


L l 


_(a) 

= X11 


X33 ' t 2 


_(a) 
— X22 



~(a) 
'X34 

~(a) 
X44 • 



(7) 



which transform as s-wave (s), o^-wave (c), <i :r 2_ J/ 2-wave (d), extended s-wave 
(e), p-wave (v x ,v y ) in the spin singlet state. Similarly, for the electron-hole 
pairs we select 

P = &11 + &22 + ^33 + <744 
P = ^11-^22 + ^33-^44 
?1 = cm - cr 33 , g 2 = cr 22 - CT 44 

m = + ^22 + ^33 + ^44 
a = ip u - (p 22 + ip 33 - <p 44 

9+ = <Pll ~ ^33 > I- = ^22 - ^44 



(8) 



They correspond to the charge density p, the charge modulation or charge 
density wave p, the ferromagnetic and antiferromagnetic spin densities to, a 

— * 

and the diagonal spin density g±. Correspondingly, we specify the source 
term for the bilinears as 



- / dr 



E 

X 



E< 



+ y^{Z 7 (x)w 7 (x) H — r 7 Z 7 (x)Z 7 (x)} — /jLw p (x)* 



..-.(2), 



+ -r p pV 2 n 2 /a 2 



with 



: (s, c, J, e,v x ,v y ,ti J 2 ) (a) , 
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The complex sources j = (j s , j c , j d , j e , j Vx , j Vy , j tl , jt 2 ) and the real sources 
I = (fx + l' p ,lp,l qi ,l g2 ,l m ,l a ,l g+ ,lg_) also depend on r and from now on we 
include the chemical potential fi in the source for p. Here V2 is the two- 
dimensional volume and we will specify the constants fp and r 7 below. 

Using the identities 

~(2)~(1) ~ ~ ~ ~ 

Xab Xcd = °caOdb + ^cb^da 

<Pab<Pcd - -<?ab&cd - 2a ad a cb (10) 
we establish the relations 

a 

Up 2 +p 2 ) + ?i + <i! = 



1/^2 ?2 

,111 —a 

2 



2 

a 

p2 + ^ +d ~(2) d ~(l) + g(2)g(l) 

+ 2(V^+V^) = 
1 f->2 -+2\ ->2 -2 sr^ 9 

-(a+m)+g + + g_ = -6^a 2 a . (11) 



The interaction term in (^) reads ^ ^ a cr 2 a and it is obvious that the decom- 
position into products of the above fermion bilinears is not unique. 

It is now straightforward to derive a partially bosonized version of the 
Hubbard model by introducing a suitable identity in the functional integral 
([]]). Inversely, the bosonized partition function 

Z = exp j-2^/3%A 2 j J D^D$*DuDu*Dw 

exp j - J dT^2(C kin + £ y + C B + £,)} 

x 

£kin = ^^(x) d T $ a {x) 
a 

- y)[ik(x, y) + ip 2 (x- 2a, y) + $4 (a, y + 2a) + ^ 4 (x, ?/)] 

+ ^2 2/) [^1 + 2a ? 2/) + 4>i( x ,y) + 4>3(Xi V + 2a) +ip 3 (x,y)] 
+ $3( x iV)$4( x + 2a ^v) +^( x ^y) + $2(x, y) (x,y -2a)] 
+ '04( x )2/)['03(^,2/) + ^3(2 - 2a, y) +ip 1 (x,y) +ip 1 (x,y- 2a)]} 
£ B = 47T 2 ^ u* (x)u[3(x) + 2tt 2 



.3 



£y = — hp(u* (x)u^\x) + U/3(x)u^\x)) — h J w 1 {x)w J (x) 
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- E (Vli^aix) + Va(x)r a (x)) (12) 
a 

can easily be transformed into a purely fermionic functional integral by per- 
forming the Gaussian integration over the complex scalar fields u and real 
scalar fields w. We choose fp = Air 2 /hp, r 7 = 4n 2 /h 2 such that the partition 
functions ([!]) and (|I2"D coincide except for the quartic interactions. Indeed, 
the four fermion interaction resulting from the bosonic functional integration 
can be more complex than in the original Hubbard model, i.e. 

Hint = ~J2 ^"fuf ~ ^7< (13) 

P 7 

Only for particular values of the real positive Yukawa couplings hp, /i 7 the 
partition function (|l"2"D is equal to the partition function (H) of the Hubbard 
modelQ namely for (cf. eq.(|TT])) 



2H S = 2H C 



H tl = Hi 



2Hd = 2H P = H„ 



3Ai 
= 6A 3 



H, 



6A 5 



qi JJ <j2 " /l 2 

H p = 3 (A 2 - A 3 ) , H p = 3 (A 2 + A 3 ) 
H 3 = 2Ai + A 2 - 3A 3 + 1 
flfli = 2Ai + A 2 + 3A 3 + 1 



3+ 



^3 



4Ai + 2A 2 + 2, 



where the parameters Aj obey 



(14) 



Ai > 
A 2 > A 3 , 
2A! + A 2 + 1 > 3A 3 . 



Vz = 1 . . . 3, 



We emphasize that the choice ( 14 ) of the Yukawa couplings is not unique 
since it depends on the three parameters Aj. Arbitrary values of Aj (within 
the allowed range) all describe the same Hubbard model. The independence 
of physical results on the values of Aj can be used as a check for the validity 



7 This holds up to an irrelevant source-independent overall normalization factor. 
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of approximations. Furthermore, a large variety of different four-fermion in- 
teractions can be described by varying the Yukawa couplings away from the 
"Hubbard values" (H). 

The symmetries R and I x as well as the translations by 2a are easily 
implemented on the space of bilinears Up,w^ and correspondingly for the 
scalar fields Up,Wj. This is not the case for the translations by a. The above 
formulation of the bosonization is therefore not optimal yet if - beyond the 
symmetries of the coarse-grained lattice - the symmetries like T x play an 
important role (as for the Hubbard model). It is easy to remedy this short- 
coming by an extension of the space of bilinears and the corresponding scalar 
fields. We introduce an additional color index for the fermion bilinears and 
the scalars by 

w 3 ^(x) = Wy(x) , w 4y (x) = T~ 1 w 7 (f) (15) 

and similar for vS a > ,w,u^ a ' Products like w^w^ are now understood as 
scalar products 

W 7 U7 7 -> - ^ W ajW aj . (16) 
a 

With these replacements^ it is straightforward to check that the partition 
function (|T2| ) with the choice of Yukawa couplings (|14]) is again exactly 
equal to the one of the Hubbard model if all sources except [i are zero. 
The translations T x ,T y are now directly implemented on the scalar fields, 
e.g. T x (w\j{x)) = W2-y(x). The same holds true for the rotation R or the 
reflection I x . 

In conclusion, we have developed an extended version of the Hubbard 
model - the colored Hubbard model - which coincides with the Hubbard 
model for special values of the Yukawa couplings (|i4[) and the sources. Par- 
ticularly simple modifications arise for x-independent and r-independent 
sources. As an example, the source 

l' 3p = h P = ~2h Ql = -\v (17) 

induces an additional energy for the occupation of the sites (m, n) with both 
m and n odd 

S„ = V dT ^ ^mn^mn- (18) 

•* m odd 

n odd 

For v —>■ oo these sites are effectively removed from the lattice and we there- 
fore deal with the Hubbard model on a non-quadratic lattice structure. 

8 The term ~ /j,w p (x) in eq. (||) becomes (/i/4) YJ a w ap (x). 
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Analytic computations for the partition function fll2|) are most easily done 
in momentum space. It is straightforward to perform a Fourier transform 
using 

/d 2 q 
exp (i{(x + z a )q + 27mTr}) i) an (q) 

exp (-i{{x + z a )q + 2-kuTt}) ^ an (q)j° 



r = 4i° , 7° = ( ; 3 ° ) ■ (19) 

Here the Matsubara frequencies are labeled by half integer n = 
±1/2, ±3/2, . . .) and the momentum integration is in the range —A < q x < 
A, —A < q y < A as appropriate for the coarse lattice with lattice distance 
2a. We choose 

a a\ _ (a a 



Zl -\2'2) ' Z2 ^ \2 ' 2 

23 = (2 ■ -2) • 24 = (-2 ■ -2) (20) 

corresponding to an expansion in the coordinates of the (m, n) lattice. This 
yields for the kinetic term 



S> 



kin J dT ^ ^ £kin 

x 

T 5Z/ -^^^an(^P^\n,q)^ bn (q) 



p(°) = -2-unT^ -2it 1 °{cos(aq x )A 1 + cos(aq y )B l } (21) 
where we use matrices (tq = I2; // = . . . 3; i, j = 1 . . . 3) 



A,= l ^ ° K B, ~( ° r " 



= 2<%£ , {A, A,} = = 2<% 

[A, -Bj] = ^ieijkBk 

B Ai = AiB = Bi , B B^ = B^B = A a . (22) 



Spontaneous symmetry breaking with "extended order parameters" like 
the antiferromagnetic spin density ~ (di— 5 2 +03—04) or d-wave superconduc- 
tivity with order parameter ~ (di ± d 2 + d 3 + <i 4 ) can be directly investigated 
in our formalism by looking for the minima of the effective scalar potential. 
The notion of the effective potential is a very powerful concept since it de- 
scribes simultaneously situations with vanishing and nonvanishing sources, 
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i.e. in addition to the Hubbard model for arbitrary p it also comprises many 
extended models. The effective potential corresponds to the effective action 
for homogeneous colored scalar fields. We define the scalar expectation values 
in the presence of sources^] 

hiZ = < Uaaix) > 







A 2 


d 




vr 2 




x) 


* i 

u a p( 




A 2 


d 




7T 2 


dJ a p( 


x) 




(x) 


A 2 


d 




7T 2 


dL a y 


(x 



\nZ = < u* a/3 (x) > 



\nZ = <w ay (x) > (23) 



with 



A 2 . A 2 

Ja/3 — -j—ja/3 , ^ay = —lay- (24) 



With the usual Legendre transform one obtains^ the effective action T 
T[u,w,ip,ip*} = -\nZ + J drJ2J2{j^J2( J li3 u af3 + J a p<p) 



x a (3 

■2 



+ X2 L ^I W ^I + Vl^a ~ Vj]a } (25) 

7 

which obeys 

dT vr 2 

— L a7 etc. (26) 



dw a y A 2 

Performing the derivatives ( f23|) in the fermionic functional integral, we can 
directly relate the scalar expectation values to the expectation values of 
fermionic bilinears 

hp ~(i) 1 j 



h 



v* - < v {2) > + — 7* 

a<i ~ 47r 2 a/3 A 2 0/3 

Wa 7 = ^ < * a7 > + ^ Ary- (27) 

In particular, if all sources except L ap (x) = A 2 p/h p vanish, the scalar w ap has 
contributions linear in the electron density n =< p > /4a 2 and the chemical 
potential 



9 The variation with respect to l ap is performed at fixed [i. 
10 We concentrate in the following on ip a — ( ip a ) = 0, ip* = ( ip* ) = 
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We will mainly be interested in homogenous expectation values and there- 
fore in the effective scalar potential which can be obtained from T for x- and 
r-independent scalar fields and vanishing fermion fields by 

U = TT/V 2 . (29) 

The ground state of the Hubbard model corresponds to the minimum of Uq 
with respect to all fields except p = \J2 a w ap given by eq. (p8|), where \i 
obeys 

hp dU , , 

* = wW m 

We are interested in possible expectation values of scalars different from p. 
Such a spontaneous symmetry breaking arises if for some range of p the 
minimum of Uq (at fixed p) occurs for a nonvanishing scalar field. 

In this paper we compute the effective potential Uq in the "mean field" 
approximation. This means that only the fermionic part of the functional 
integral (|12"D is performed in a homogenous "background" u = u, w = w. 
This integral is Gaussian, and we can write the mean field expression for Uq 
as 

Uq = U cl + AU 

U d = A 2 Yl u *aP u ^ + Y S W ^ W ^ + "XT V 2 

P a 7 a P 

72. 



AU 



\ T Hj ^IndetP^g). (31) 



Here P(q) is a 16 x 16 matrix (including spinor indices) for the inverse fermion 
propagation in the presence of scalar background fields. It is defined by the 
part of the action quadratic in the fermion fields 



d 2 q 7T 




iPt n (-q)P(n,q)MQ) (32) 

(33) 



where the second term reflects the influence of the background through the 
Yukawa couplings. We explore here the dependence of the effective potential 
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on the charge density p, d-w&ve pair condensation d and antiferromagnetic 
order parameter a. We therefore take 



wi p = w 2p = w 3p = w Ap = p 
uid = u 2d = u 3d = u Ad = d 

wi a = -w 2a = w 3a = -w 4a = a (35) 



and find0 



P = -ihpf ° -^\-iha( "-WW 

d*[cos(aq x )Ai — cos(aq y )Bi] 
v d , y°[cos(aq x )Ai — cos(aq y )B 1 ]'y o1 

P T . (36) 



With G = diag(l, -27° ) and P = GPG T one obtains 

In det P = In det P = - In det (PB P T B ) 

= In detg{(27mT) 2 + [2t cos(aq x )Ai + 2tcos(aq y )Bi + h p p] 2 
+ h 2 a ad + 2h p h a paA 3 ® f 

+ h d d*d[cos(aq x )Ai — cos(aq y )B 1 } 2 }. (37) 

It is easy to see that AU depends only on the invariants 5 = d*d, a = aa 
and p. Up to an additive (T-dependent) constant one finds 



AU = --T > / T^rrtrg In <i 1 



2 (27T) 



d 2 g 



\2 



+ f [2tcos(ag x .)Ai + 2tcos(aq y )Bi + /i p p + /^v 7 "^ ® r 3 ] 

+ fc^ [cosfag^Ai - cos(ag 2/ ) J B 1 ] 2 ) /(2vmT) 2 } (38) 



C/ d = 2A 2 a + 4A 2 5 + 2AV + ^^, 

and, evaluating the Matsubara sum and the trace, finally 

a, ,2 r ai 



P V 7 £i ,£j6{-l,l} 

In cosh ^7^^/ (V> + £i\pti 2 ~{cx + £jC y ) 2 + h 2 a a^j + h 2 5(c x - £jC y ) 2 ^j 

with c x = cos(aq x ), c y = cos(aq y ). 

11 The choice of 7 in (|l9|) is not crucial for this calculation - any orthogonal 7 will do. 
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For large temperature the fluctuation contribution AU is suppressed ~ 
T _1 . The minimum of Uq therefore occurs for all p at a = 0, 5 = 0. As T is 
lowered, the fluctuations tend to destabilize the "symmetric minimum". In 
particular, the fluctuation contribution to the mass term for a is negative for 
not too large p 2 and the one for d is negative for all p 

AM 2 a = 2^-AU ]a=s=0 

= -2hlTj2 J ^^ 2 -2hyA 3 p- 2 A 3 P^} 



d 



AM? = ^-AUi 



' h d T Yl / ^tr 4 {(cos(ag :c )A 1 -cos(ag 2/ ) J B 1 ) 2 p- 2 } (39) 



with 



P 2 = (2vmT) 2 + (2t cos(aq x )A 1 + 2t cos(aq y )B 1 + h p pf (40) 

and tr 4 the trace in color space only. These contributions should be compared 
with (Mi 0) ) 2 = (Mf ] ) 2 = 4A 2 . The zeroes of P 2 for T = correspond to 
the Fermi surface (|3]) with shifted chemical potential /z e fr = h p p. (Neglecting 
contributions from AU eqs. (|30| ) and (|3T1) imply p cS = p.) We recall that 
the momenta are restricted to the range corresponding to the coarse lattice 
\Qx,y\ < 7r/(2a) = A. On the other hand we have now possible zeroes for 
different linear color combinations. Noting that the eigenvalues of A\ and B\ 
are ±1, they precisely correspond to the original Fermi surface - the original 
zeros in the four ranges \q x , y \ < A, A < \q x , y \ < 2A, \q x \ < A, A < \q y \ < 2A 
and A < \q x \ < 2A, \q y \ < A appear now for different color combinations 
in the range \q x ,y\ < A. Due to these zeros one finds linrr^oAMj — > — oo 
for not too large p and similar for AM 2 in the appropriate range of p. This 
clearly indicates spontaneous symmetry breaking with d-w&ve superconduc- 
tivity or/and antiferromagnetic order parameter. Note that in contrast to its 
derivatives the potential is not singular for T —>■ 0. Since for large a and 5 
Uq grows ~ (a, 5), the minimum occurs necessarily for finite a ^ or 5 ^ 
if M 2 = 4A 2 + AM 2 or Mj = 4A 2 + AMj become negative. 

The spontaneous symmetry breaking cuts off the singularity near the 
Fermi surface or reduces its strength. An antiferromagnetic expectation value 
typically produces a gap for the fermionic fluctuations. For a > 0, 5 = 0, p = 
this can be seen from a search for possible zeroes of dets in eq. (|37|) for 
T = 0. On the other hand, for a = 0, 5 > the condition dets = requires 
cos(ag x ) = ±cos(aq y ) = ±h p p/(4t). In the superconducting phase the sin- 
gularity therefore only occurs for special points in momentum space instead 
of a whole Fermi surface. As a consequence, the momentum integrations for 
the bosonic mass terms (similar to eq. (|39"D) remain finite even for T — » 0. 
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0.25 




p/Vu 

Figure 2: The T — p phase diagram for h p = h a = hs = VlOU with symmetric 
(SYM), antiferromagnetic (AF) and superconducting phase (SC). In the re- 
gion marked by the bold line the phase transition into the antiferromagnetic 
phase is of first order; all other phase transitions are of second order. 



For vanishing sources j a , jd the minimum of U obeys the "field equations" 

< 41 > 

tr 8 {( h pP a- 1/2 A 3 <g> r 3 + h a )P~\a, 5)}) = 

i^-\^Sw? (42) 

tr 8 {[cos(a^) - cos(aq y )Bo] 2 P p 2 (a,5)}) = 



dU 2 M o ^ oa2 [ d 2 q 

da da 



dU dU j ax2 f d 2 q 

dd* 85 



with 



P 2 Ja,5) = P 2 + 2h p h a pVaA 3 ® r 3 



+ h 2 a a + h 2 d 8 [cos(aq x ) — cos(aq y )B ] 2 . (43) 

One always has the symmetric solution a = 0, 5 = which corresponds to 
a local minimum if M a > 0, M d > and to a maximum or saddlepoint 
otherwise. Consider next solutions with a = 0, 5 ^ which require 



/d 2 q { 4 
^-^tr 4 | [cos(a^) - cos(aq y )B ] 4 P f 



4 o-2 
P 
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0.1 0.2 0.3 0.4 0.5 



p/Vu 

Figure 3: The T — p phase diagram for h p = h a = h$ = 2y/10U with sym- 
metric (SYM), antiferromagnetic (AF) and superconducting phase (SC). In 
the region marked by the bold line the phase transition into the antiferro- 
magnetic phase is of first order; all other phase transitions are of second 
order. 

" ( p p + h d$ [cos(aq x ) - cos(aq y )B ) 2 ) 1 J 

= 4A 2 + AM ( 2 = M d . (44) 

2 

Solutions with 5 > exist only for M d < and 5 vanishes as the mass term 

2 

M d approaches zero from below. One concludes that the transition from the 
symmetric phase (a — 0, 5 = 0) to a possible superconducting phase without 
antiferromagnetism (a = 0, 5 > 0) is of second order. 

We have analyzed the phase diagram for different Yukawa couplings 
numerically. Due to the free parameters Aj in eq. ( |I~4T) the Yukawa cou- 
plings are largely undetermined. They only have to obey the inequalities 
h 2 p > 0, h\ > 0, h 2 a > 0, h 2 a > tt 2 U/3 + h 2 p /3 - 2h 2 d /3. For example, 
\ 1 = A 3 = 1/2, A 2 = 1 leads to h 2 d = h 2 = h 2 a = n 2 U/2. Because 
of our meanfield approximation, the partition function becomes dependent 
on the particular choice of the parameters Aj. We investigate the cases 
h p = h a = h 5 = VlOU (fig. |), h p = h a = h s = 2 v / 10f7 (fig. 1), K = VlOU, 
h p = h s = 2 V / T0Z7 (fig. |), and h 5 = VlOU, h p = h a = 2 v / 10f7 (fig. |). We 
choose t/U = 1 and investigate the phase diagram in the p/VU-T/U -plane. 
Expressed in the variables t/U, T/U, pj VU, our results do not depend on U 
and the lattice distance a, as discussed in the beginning. As we increase all 
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Figure 4: The T — p phase diagram for h a = Vl~0Z7, h p = hs = 2\A0£7 with 
symmetric (SYM) and superconducting phase (SC). The phase transitions is 
of second order. 




SYM 



0.1 0.2 0.3 0.4 0.5 

p/Vu 

Figure 5: The T — p phase diagram for h 5 = y/lQU, h p = h a = 2v / 10C7 with 
symmetric (SYM) and antiferromagnetic phase (AF). In the region marked 
by the bold line the phase transition into the antiferromagnetic phase is of 
first order, otherwise of second order. 
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three Yukawa couplings simultaneously, the antiferromagnetic phase domi- 
nates over the superconducting phase (compare figs. ^andQ). An interesting 
result of the mean field analysis is the appearance of a phase transition of first 
order into the antiferromagnetic phase for small T/U and for high values of 
pj \fu . The phase transition between the symmetric and the superconducting 
phase remains of second order. Both results were anticipated when examining 
the above formulae analytically. If we increase hd/VU compared to h a / y/U 
the superconductivity phase dominates for low T/U, whereas in the opposite 
case it is the antiferromagnetic phase. This is illustrated in figures f| and [5[ 

We note that for negative t our results apply if the antiferromagnetic 
condensate a is replaced by the ferromagnetic condensate m. Furthermore, 
small disturbances can easily be taken into account by source terms. For 
example, an interaction between spin and angular momentum will explicitely 
break the continuous SU(2) invariance and typically amount to a source term 

la Or lifi- 

In conclusion, the mean field approximation for the colored Hubbard 
model can give a qualitatively reasonable picture of the phases in high T c 
superconductors. On the other hand, the shortcomings of this approxima- 
tion are also apparent from the figures. All phase diagrams in figs. 0, |] 
and |5| correspond to different mean field approximations for the same model. 
It is impossible to resolve this ambiguity within the mean field approxima- 
tion without additional input on the selection of the Yukawa couplings. The 
reason is the neglect of fluctuations of the bosonic fields. Only if these are 
included, the different equivalent choices of the Yukawa couplings should lead 
to the same physical results. The differences between the figures reveal the 
importance of the neglected bosonic fluctuations, at least for some choices of 
the Yukawa couplings^. 

The inclusion of the bosonic fluctuations is a complex problem which can 
be attacked by means of nonperturbative renormalization group equations |7j . 
Studies for similar QCD-motivated models of fermions with Yukawa coupling 
to scalars have already been carried out successfully || . One of the dominant 
effects will be the scale dependence of the Yukawa couplings. It is conceivable 
that this running is dominated by partial infrared fixed points for ratios of 
Yukawa couplings. For large couplings, as relevant here, such partial fixed 
points would be approached fast. In this case the "memory" of the initial 
choice of Yukawa couplings could be erased rapidly and unambiguous physical 
predictions become possible. 

A second important ingredient is the appearance of Goldstone bosons for 
(a) 7^ or (d) ^ 0, corresponding to flat directions in the effective poten- 
tial (pTtp . For a superconducting condensate (d) the ?7(l)-symmetry would 

12 It is conceivable that an "optimal choice" of the Yukawa couplings minimizes the 
impact of the bosonic fluctuations. 
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be spontaneously broken and the question arises if this is self consistent. For 
a large correlation length £, i.e. £T ^> 1, one expects that the dominant 
fluctuations near a second order phase transition are well described by an 
effective dimensional reduction to two dimensional classical statistics. The 
Mermin-Wagner theorem then suggests that the Goldstone boson fluctua- 
tions prevent a continuous symmetry from being spontaneously broken. In 
the case of a [/(l)-symmetry the natural solution to this puzzle is a second 
order phase transition of the Kosterlitz-Touless type: only a renormalized ex- 
pectation value differs from zero, whereas the wave function renormalization 
will lead to a vanishing expectation value for the unrenormalized scalar field 
H . This reconciles the Mermin-Wagner theorem with the existence of Gold- 
stone bosons and superconductivity in presence of electromagnetic fields. 

For a possible "antiferromagnetic phase" the nonabelian interactions be- 
tween the Goldstone bosons of the effective two dimensional model have a 
tendency to push the minimum of Uq towards a = and to make dUo/da 
positive 0. If only the nonabelian Goldstone bosons are present in the ef- 
fective long distance model their fluctuations would destroy the nontrivial 
minimum of the potential. One may therefore speculate about a new type 
of low temperature phase, which is characterized by the presence of mass- 
less Goldstone bosons as well as massless fermions. Alternatively no true 
antiferromagnetic phase with Goldstone bosons may occur. For all practical 
purposes the physics nevertheless will look qualitatively similar to the phase 
transition in the mean field approximation: the effects from Goldstone fluc- 
tuations are only logarithmic in ratios of mass scales and would be cut off by 
a small S{7(2)-breaking disturbance inducing a mass term for them. Simple 
scale considerations suggest that the first order transitions to the antifer- 
romagnetic phase are probably not affected substantially by the Goldstone 
fluctuations, except for the endpoints. Particularly interesting is the triple 
point in fig. ^ where the three phases meet. By continuity of the second order 
lines one expects five massless scalar excitations at this point. 

We emphasize that quite generally the possible second order phase tran- 
sitions between the symmetric and some other phase belong to new interest- 
ing universality classes. Long range fermion fluctuations without a gap are 
present in the symmetric phase and therefore also at the phase transition. 
They influence the critical exponents and other universal properties. 

We hope that our formulation of the colored Hubbard model will be a 
good starting point for a quantitative renormalization group study of all these 
interesting questions. 
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